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1 Introduction 



The central concept of the random matrix theory as envisioned by E. Wigner is the general hy- 
pothesis that the distributions of eigenvalue spacings of large complicated quantum systems are 
universal in the sense that they depend only on the symmetry classes of the physical systems but 
not on other detailed structures. The simplest case for this hypothesis is for ensembles of large but 
finite dimensional matrices. The general hypothesis in this setting thus asserts that the eigenvalue 
spacing distributions of random matrices should be independent of the probability distribution of 
the ensemble, up to scaling. This is generally referred to as the universality of random matrices. 
In this paper we will focus only on the bulk behavior i.e., on eigenvalue distribution in the interior 
of the spectrum, although similar questions regarding the edge distribution are also important. 

Over the past two decades, spectacular progress (see, e.g., [5, 10, 11, 24, 25, 9, 22] and [2, 8, 9] 
for a review) on bulk universality was made for classical invariant ensembles, i.e., matrix models 
with probability measure given by jZ where is the size of the matrix 1/ is a 

real valued potential and Z is the normalization. It is well-known that the probability distribution 
of the ordered eigenvalues of H on the simplex determined by Ai < • • • < Ajv is given by 

^ = MA, -A.), (1.1) 

k=\ l<i<j<N 

where the parameter /3 = 1,2,4 is determined by the symmetry type of the matrix, corresponding 
respectively to the classical orthogonal, unitary or symplectic ensemble. With P taking these special 
values, the correlation functions can be explicitly expressed in terms of polynomials orthogonal 
to the measure e~'^^^^'/^. Thus the analysis of the correlation functions relics heavily on the 
asymptotic properties of the corresponding orthogonal polynomials. In the pioneering work of 
Gaudin, Mehta and Dyson (see [23] for a review), the potential V is the quadratic polynomial 
V{x) = and the orthogonal polynomials are the Hermite polynomials for which asymptotic 
properties are well-known. The major input of the recent work is the asymptotic analysis of the 
orthogonal polynomials w.r.t. the measure e"'^^^^-'/^ for general classes of potentials. The formulas 
for orthogonal and symplectic cases, i.e., /? = 1,4, are much more difficult to use than the one 
for the unitary case. While universality for /? = 2 was proved for very general potential, the best 
results for /3 = 1,4 [9, 21, 26] are still restricted to analytic Y with additional conditions. 

For non-classical values of /3, i.e., /3 ^ {1,2,4}, one can still consider the measure (5.47), 
but there is no simple expression of the correlation functions in terms of orthogonal polynomials. 
Furthermore, the measure (5.47) does not arise from mean-field type matrix models like Wigner 
matrices with independent entries. Nevertheless, /i is a Gibbs measure of particles in R with a 
logarithmic interaction, where the parameter /3 is interpreted as the inverse temperature and a 
priori can be an arbitrary positive number. These measures are called general /3-ensembles. We 
will often refer to the variables \j as particles or points and the system is called log-gas. It was 
proved [12] that in the Gaussian case, i.e., when V is quadratic, the measure (5.47) describes 
eigenvalues of tri-diagonal matrices. This observation allowed one to establish detailed properties, 
including the local spacing distributions of the Gaussian /3-ensembles [27]. 

Gibbs measures in the continuum with long range or singular interactions are notoriously hard 
to analyze since they are very far from the perturbative regime. For non-classical values of /3, 
and if we are not in the Gaussian case ViX) = A^, no simple explicit formula is known to express 
the correlation functions in terms of orthogonal polynomials, and one cannot rely on any explicit 
known matrix model. In this paper we undertake the direct analysis of the Gibbs measure and 
we prove the universality for invariant models for any /? > 0. In other words, we will prove that 
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the local spacing distributions of (5.47) are independent of the potential V for certain class of V. 
There are two major ingredients in our new approach. 

Step 1. Uniqueness of local Gibbs measures with logarithmic interactions. The main result in this 
step asserts that if the particles are not too far from their classical locations then the spacing 
distributions are given by the corresponding Gaussian ones (We will take the uniqueness of the 
spacing distributions as our definition of the uniqueness of Gibbs state). More precisely, denote 
by p the limiting density of the particles under the measure /j,*^^^ (5-47) as TV — >■ oo. Let 7j = ^j^N 
denote the location of the j-th point under p, i.e., 7^ is defined by 



We will call 7^ the classical location of the j-th particle. The basic assumption is the following: 
Assumption A. For some h < and any a > 0, there exists £0 > such that 



for large enough N and any k e [aN, {1 ~ a)N]. 

Under this assumption (under some minor and easily verifiable assumptions near the edges of the 
limiting measure) , we will prove that the spacing distributions of are given by the corresponding 
Gaussian model with V{x) = . We will use the Gaussian case as our reference ensemble only for 
the convenience of definiteness. In fact, no detailed properties of the Gaussian measures are used 
in the proof and any other reference ensemble would have worked as well. Furthermore, in this 
step we make no assumption on the convexity of V , which is needed in the next step. 

Step 2. Particle location estimate. The second step is to verify Assumption A. For non-classical 
/3, Assumption A is only proved for b near one [20, 25, 21] for analytic potential V under certain 
constraint. This is far from sufficient to complete Step 1. We will prove Assumption A for all /3 > 
under the assumption that V is convex and analytic. Our method uses the following three ideas: 
(1) The analysis of the loop equation in [20, 21, 26] to control the density. (2) The logarithmic 
Sobolev inequality guaranteed by the convexity of V . (3) A multiscale analysis of the probability 
measures of invariant ensembles. We note that the assumption of analyticity on V is needed only 
for using the loop equation in (1). 

The basic idea of our proof is to use the following tool from [14]: For two probability measures 
/i and uj define the Dirichlet form by 



Then the difference of the local spacing distributions of the two measures is negligible provided 
that the Dirichlet form per particle is sufficiently small in the large N limit [14]. Notice that if we 
used the relative entropy of the two measures, then the uniqueness of the Gibbs measures would 
require the total entropy, which is an extensive quantity, to be small. To apply this Dirichlet form 
inequality, we first localize the measure by fixing \j for j outside, say, the interval [L + 1,L + K] 
for L in the bulk and K = for some A: > 0. We will call these data of \j outside the interval 
[L + 1,L + K] the boundary condition. We then compare this measure to a local Gaussian /?- 
ensemble with a fixed boundary condition by showing that the Dirichlet form per particle of these 
two measures is small for typical boundary conditions w.r.t. /x. 




(1.2) 



P^(«)(|Afe-7fc| < > l-exp(-iV^'') 



(1.3) 
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Our approach shares some philosophy from the recent method on the universahty of Wigner 
matrices [14, 15]. In this approach, the key condition to estabhsh is 

Assumption III. There exists an o > such that wc have 

^,..j^Y.^^,-li?<CN-^-^<^ (1.4) 
i=i 

with a constant C uniformly in N . Here is the law given by the Wigner ensemble. 

Under this assumption, a strong estimate on the local ergodicity of Dyson Brownian motion 
(DBM) was established in [14, 15]. DEM [13] establishes a dynamical interpolation between Wigner 
matrices and the invariant equilibrium measure ^. This estimate then implies the universality of 
Wigner matrices. Thus the main task in proving the universality of Wigner matrices is reduced to 
verifying Assumption III. 

There are several similarities between the method used for the universality of Wigner matrices 
[14, 18] and the current proof for /3-cnscmbles: (i) Both rely on crude estimates such as (1.3) and 
(1.4) on the location of the eigenvalues to establish the local spacing distributions are the same as 
in the Gaussian cases, (ii) Both use estimates on the Dirichlet form to identify the local spacing 
distributions, (iii) The main model dependent argument is to prove these crude bounds on the 
eigenvalues. The precision of these a-priori estimates on the eigenvalues is weaker than the local 
spacing, but better than previously known results on eigenvalue locations: we have to develop new 
methods to prove (1.3) and (1.4). 

There are, however, substantial differences between the proofs of universalities for Wigner 
and /3-ensembles. First, since the /3-ensembles are already in equilibrium, there is no dynamical 
relaxation mechanism to exploit and the local statistics need to be identified directly without 
dynamical argument. Second, we obtain the crude estimate (1.3) by a method completely different 
from the Wigner matrices, as there is no underlying matrix ensemble with independent entries to 
analyze. The accuracy result we obtain by this new method is actually optimal, i.e. (1.3) will be 
shown to hold for any b > 0. 



2 Statement of the main result 

Consider a probability measure 

^,f^ = ^i^)[AX) = ^ n |A.-A,|^^^"''^''^'^'dAl...dA^., (2.1) 

^ l<i<j<Af fc=l 

where A = (Ai, . . . , Aat), Ai < • • • < Xn- Here the inverse temperature satisfies /3 > and the 
external potential V is any convex real analytic function in R, and such that 

^7 = f inf r'(x) > 0. (2.2) 

For such a convex potential, as noted in the next section the equilibrium measure, denoted by 
p(s)ds, is supported on a single interval [A, B]. In the following, we omit the superscript N and 
we will write /i for We will use and to denote the probability and the expectation 

with respect to ^. 
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The Gaussian case corresponds to ^(A) = A^; the expectation with respect to this Gaussian 
measure will be denoted by Ecauss , and the equilibrium measure is known to be 



Psc{E) — v/(4-i?2)+, 

the semicircle density. The Gaussian case includes the classical GUE, GOE and GSE ensembles 
for the special choice of /3 = 1,2,4, but our result holds for all /3 > 0. 

Now we state our main theorem which will be proven at the end of Section 4: 

Theorem 2.1 Assume V is any real analytic function with inf V"{x) > 0. Let j3 > 0. Consider 
the /3-ensemble jjL = fipy. Let G : R — > R &e a smooth, compactly supported function. Let 
E e (A, B) lie in the interior of the support of p, and similarly let E' € (—2, 2) he inside the 
support of psc ■ Define L and L' by 



L 

N 



p{x)dx, 
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E' 



Psc{x)dx. 



Fix a parameter K ~ N'^ where < fc < 5 is an arbitrary constant. Let L and I' he two intervals 
of natural numbers, I = [L + 1, L + A'] , /' = [L' + 1, A' + AT] with length A' = |/|. Then 



lim 

TV— s-oo 



Kp{E) 



■lei 



N{\ - A,+i) 
P{E) 



Psc{E') 



0, 



(2.3) 



i.e. the appropriately normalized particles gap distribution of the measure fi/3_v o-t the level E in 
the bulk of the limiting density asymptotically coincides with that for the Gaussian case and it is 
independent of the value of E in the bulk. In particular the gap distribution is universal. 

Remark. The same result (with the same proof) holds for higher order correlation functions of 
particles gaps. More precisely, fix n > 1 and an array of positive integers, m = {mi, TO2, ■ • ■ , fn„) e 
N^L. Let G : M" — !• R be a bounded smooth function with compact support and wc define 



Gi,inW ■ = 



1 



p{Ey 



-G 



iV(A 



i-\-mn 



P{E) 



P{E) 



P{E) 



Then, under the conditions of Theorem 2.1 and using its notations, we have 



lim 



E„ 



E: 



'Gauss 



0, 



(2.4) 



(2.5) 



iei iei' 

where Gi_m is defined exactly as Gi,m but p{E) is replaced with psc{E'). 

The limit (2.5) can be reformulated as the convergence of the correlation functions. Let pi^^ 
denote the n-point correlation function of the measure p. = p'j^y defined by 



(xi 



fl{x)dxn+i . . .dxN, (2.6) 
where p is the symmetrized version of p given in (2.1) but defined on R^ instead of the simplex: 

AW(dA) = ^MdA('^)), 
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where A^'^^ = (A^(i), . . . , A^(7v)), with A<^(i) < • • • < A<^(jv)- 

From (2.5) we have the convergence of the correlation functions, stated as the following corol- 
lary. Since the proof is a standard argument and it is essentially identical to the one given in 
Section 7 of [15], we omit it. 

Corollary 2.2 Under the assumption of Theorem 2.1 and with the same notations, for any smooth 
test functions O with compact support and for any < fc < we have, with s := N^^^^ , that 



J- 



lim / dai • • • dan 0(ai, 



, an) 



E+s 



dx 1 



1 



^+^dx 



(N) 



■Pa 



E-s 2sq{E] 
oil 



NQsc{E) ' 



Nq{E) '■ 



NQsc{E) 



an 



Nq{E)) 
0. 



The local statistics of the X'^s in the Gaussian case have been explicitly computed by Gaudin, 
Mehta and Dyson (see, e.g., [23]) for the classical value /3 S {1, 2, 4}. For general /? > 0, there is an 
explicit description in terms of some stochastic differential equations, the Brownian carousel [27]. 

Theorem 2.1 will be proved in two steps as explained in the introduction. For logical reasons, 
we will first present Step 2 on particle location estimates in Section 3 and then Step 1 on the 
uniqueness of Gibbs measure in a finite interval in Sections 4 and 5. 



3 Optimal accuracy for particle locations 

Along this section, we assume that V satisfies the same conditions as in Theorem 2.1. Let the 
typical position 7^ be defined by 

/7fc ^ 
-00 ^ * 

Moreover, all constants in this section depend on the potential V, which is fixed. In the following, 
we wiU denote |x, y] = N n [a;, y], 

The purpose of this section is to prove that accuracy holds for the measure fi at the optimal 
scale 1/A^, in the following sense. 

Theorem 3.1 Take any a > and e > 0. There are constants 8,C\,C2 > such that for any 
N >l andke faN, (1 - a)Nl 

After some initial estimates relying on large deviations results, the proof consists in comparing 
fj, to some locally constrained measures for which better concentration estimates can be proved 
for the differences between particles. This measure is related to the pseudo-equilibrium measure 
in [14], but has distinctly different properties. Iterations of these comparisons will give optimal 
accuracy. 



3.1 Initial estimates 

The purpose of this paragraph it to prove the following crude estimate. It will be the initial step 
in the induction of Subsection 3.3. 
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Proposition 3.2 For any a,s > there are constants ci,C2,S > such that for any N and 
k e laN, (1 - a)Nl 

(|Afe - 7fc| > N-i+') < cie-^^^'. (3.1) 

This result is a direct consequence of tire following equation (3.12) and Corollary 3.5, whose 
proofs arc the purpose of this section. We first state well-known facts about the equilibrium 
measure. 

For convex analytic potential V satisfying the asymptotic growth condition (2.2) (or even with 
weaker hypotheses on V, see e.g. [6, 1]), the equilibrium measure p{s)ds associated with {fJ.^^^)N>o 
can be defined as the unique minimizer (in the set of probability measures on M endowed with the 
weak topology) of the functional 

/(z/)= J V{t)dv{t)- jj \og\t - s\du{s)dv{t) 

if / V{t)dv{t) < oo, and = oo otherwise. Moreover, p has the following properties: 

(a) The support of p is a single interval [A, B] . 

(b) This equilibrium measure satisfies 

ino - J (3,2) 

for any t e {A, B). 

(c) For any t G [A,B], 

p[t) dt = K{t)^ {t - A){B - t)\A,B] dt, (3.3) 
where r can be extended into an analytic function in C satisfying 

27tJa z-t ^{t-A){B-t) 
In particular, for convex V , r has no zero in R. 

It is known that the particles locations cannot be far from its classical location [4, 26]: for any 
e > there are positive constants C, c, such that, for all > 1, 

{3k e ll,Nl I |Afc - 7fc| > £) < Ce-^^\ (3.5) 

In order to have density strictly in a compact support, for given R > 0, define the following variant 
of /i^^^ conditioned to have all particles in [— -R, R]: 



1 ^ 
^(^■«)(dA) = -^ n |A.-A,|^n^"''^''''^^l|A.I<«dAi...dA^. (3.6) 



l<i<j<N k=l 



Let p\f^'^'^ denote the marginals of the measure p'^^'^\ i.e. the same definition as (2.6), but with 
replaced by p'-^'^\ 
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Then Lemma 1 in [6] states that under condition (2.2) there exist some R > and c > 0, 
depending only on V, such that for any |xi |, . . . , < R 

p[!^''^\xi,...,Xk) - pi^\xi,...,Xk) < pf'^\xi,...,Xk)e~''^, (3.7) 

and for |xi|, . . . , |a;j| > R, \xj+i\, . . . , \xk\ < R, 

(.xi,...,Xfc)<e-^^t.iog|.,| (38) 

The last type of estimates we need are concentration and accuracy of the particles location at scale 
N~^^^, in the bulk. Concentration is a simple consequence of the Bakry-Emery convexity criterion 
for the logarithmic Sobolev inequality ([3], see also [2]): define "H by ^(dA) = -^e~^^^^^dX, and 
assume 

\'^H>aldN (3.9) 

in the sense of partial order for positive definite operators. Then p, satisfies a logarithmic Sobolev 
inequality with constant 2 / (crN) : for any probability density / we have 

EM/log/< J^E^|Vv/7|2. (3.10) 

It is well-known that the logarithmic Sobolev inequality implies the spectral gap and, together 
with Herbst's lemma, it also implies that for any fc € |1, NJ and a; > 

P^(|Afc-E^(Afe)| >.t) <2e-'^^-'/^ 
In our case where p is defined by (2.1), for any v G 

i<j i 

where vj is defined in (2.2). Tlius there is a constant c > such that for any k 

(|Afc - E^(Afc)| >x)< 2e-^^"', (3.12) 

i.e. concentration at scale 1/\/N holds. We now prove that accuracy at the same scale holds inside 
the bulk. 

The proof of the following lemma is based on an argument in [20] for the polynomial case. In 
the form presented here, it follows very closely the proof in [26] for the analytic case except that we 
use the logarithmic Sobolev inequality to have a more precise estimate. We now introduce some 
notations needed in the proof. 

• toat is the Stieljes transform of pj^'' (s)(ds), evaluated at some z with Im(z) > 0, and m its 
limit: 

mA.(z)=Ejlf^^ ) = / ^pW(i)dt, m(z)= / ^p(t)dt. 

It is well-known that uniformly in any {Im(2:) > e}, e > 0, [toat — m| (see e.g. [2]). 
Along the proof of the next Lemma 3.3 we will see that this convergence holds at speed 1/A''. 
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• s{z)^ -2r{z)^/{A- z){B - z), where the square root is defined such that 



f{z) = \/{A- z){B - z) ~ z as z o 
hniz) is an analytic function defined by 

JR Z — I 

finally, cn{z) = -^kN{z) + ;^ (| - l) rn'j^i.z), where 

fcAr(z) = var^ (XIt^X ) ■ 




Here the var of a complex random variable denotes var(X) = E(X^) — E(X)^, i.e. without 
absolute value unlike for the usual variance. Note that |var(X)| < E(|X — E(X)p). 

The equation used by Johansson (which can be obtained by a change of variables in (2.1) [20] or 
by integration by parts [26]), is a variation of the loop equation (see, e.g., [19]) used in physics 
literatures and it takes the form 

(mjv — TiiY + s{mN — rn) + bjy ^ cn- (3.13) 

Equation (3.13) expresses the difference rajq — m in terms of (mjv — m)^, bjq and cn- In the 
regime where [tojv — '^^| is small, we can neglect the quadratic term. The term is the same order 
as Im^r — m| and is difficult to treat. As observed in [1, 26], for analytic this term vanishes 
when we perform a contour integration. So we have roughly the relation 

(3.14) 

where we dropped the less important error involving m'j^{z)/N due to the extra 1/A^ factor. In 
the convex setting, the variance can be estimated by the logarithmic Sobolcv inequality and we 
immediately obtain an estimate on to^v ~ ™- We then follow the method in [16] to use the Helffer- 
Sjostrand functional calculus to have an estimate on the particle locations. Although it is tempting 
to use this new accuracy information on the particle locations to estimate the variance again in 
(3.14), this naive bootstrap is difficult to implement. The main reason is, roughly speaking, that 
the particle location estimate obtained from knowing only the size of tojv — m is not strong enough 
in the bootstrap. The key idea in this section is the observation that accuracy information on 
particle locations can be used to improve the local convexity of the measure fj, in the direction 
involving the differences of particle locations, see Lemma 3.8. Now we are able to complete the 
bootstrap argument and obtain a more accurate estimate on — m. Since this argument can be 
repeated, we can estimate the locations of particles up to the optimal scale in the bulk. 

Lemma 3.3 Let 6 > 0. For z = E + b-j with A + d<E<B — S assume that 

l^kN{z)^0 (3.15) 
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as N oo uniformly in t] > N for some < a < 1. Then there are constants c,e > such 
that for any N-'^+'' < i] < e, A + S < E < B - S, 

|tojv(z) < c + ^fcAr(z)^ . (3.16) 

Proof. First, for technical contour integration reasons, it will be easier to consider the measure 
(3.6) instead of ^^^^ here. More precisely, define 



Then it is a direct consequence of (3.7) and (3.8) that there arc constants c > and R > such 
that uniformly on ry > N~^^ (or any power of N), 

|m^)-mAr| = 0(e-^^), l^^^ - fcAr| = 0(e"^^). (3.17) 

Consider the rectangle with vertices 2R + i7V-i°, -2R + iTV^io, -2i? - iiV^i", 2R - iN-'^°, call £ 
the corresponding clockwise closed contour and C the one consisting only in the horizontal pieces, 
with the same orientation. From (3.13), we obviously have, for z ^ C , 

1 f (m^(e) - mjOf + sjOimNiO ~ m(g )) + 6^(0 - cn{0 _ ^ 



Note that the above expression makes sense for large enough N, because then r has no zero on C. 
Using (3.17), this implies, for 77 > N~^, 



Now, as Pi ' and p are supported on [— i?, i?], m)^ — m and cj^ are uniformly 0(1) in the 
vertical segments of C. Consequently, from the above equation 

As and r arc analytic inside C, for z outside £ we get 



^ /• (^^(0 - MO)' + -KOl^^riO - ^(0) - c'n 'iO ^ ^,^-10^ 
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Remember we define f{z) = y^(A — z){B — z) uniquely by f{z) ^ z as z — > oo. Moreover, 



'AT 



77i|(z) = 0{z ^) as |z| oo because p and p^^'^'' are compactly supported 



\m)^'{z) - 771(2)1 = 



z-t 



dt 



(t)) 



o 



dt 



O (Z-) 



Consequently, the function s{mj^'' — m^/r = ~2f{m]y ' — to) is 0(z ^) as |z| —J- 00. Moreover, it 
is analytic outside C, so the Cauchy integral formula yields 



de=-2/(z)(TO^)-m)(z), 



proving 



-2f{z){m\^\z)~m(z))^- — 



1 f i^'-NHo-mior-c'j^'io 



d^ + OiN~'"). (3.18) 



Consider now the following rectangular contours, defined by their vertices: 

Ci : R + e + ie, —R — e + ie, — i? — e — ie, i? + e — ie, 

C2: R + 2e + 2i£, -i? - 2e + 2i£, -i? - 2e - 2ie, i? + 2£ - 2ie, 



(3.19) 



where £ > is fixed, small enough such that all zeros of r are strictly outside £2- For z inside £2 
and Im(z) > iV"-'^, by the Cauchy formula, equation (3.18) implies that 



2s(z)(miv''^^(z)-TO(z)) 
= -(to(^^)(z)-to(z))2 + 4^)(z)-^ 



(^!v-^(e)- »;(e))- - 4-^(0,^^0(^-10). (3 ,0) 



In the above expression, if now z is on Ci, |z — ^| > e, and on £2 |^| is separated away from zero 
by a positive universal constant. Moreover, c^j^\^) can be bounded in the following way: for any 



constants ai, . . . , & [—R — e, R + e], 



1 

iV2 



\k=l ^ 



^ 1 



1 



because for any /c. we have |^ — Xk\ > £, \£, — ak\ > £■ Now, choose ak = Ep(AA;). By (3.8), for 
large enough N any ak, 1 < k < N, is in [—R — e,R + e] indeed. Moreover, by (3.7), 



|E^(«.H) (|Afc - afcl' ) - (|Afc - ak\^) \ < e-^^E^(«,H, {\Xk ~ akf) , 
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and, by the spectral gap inequality for /i, (|Afc — afcp) = 0{N ^). This together proves that 

^N^^O is 0(A^~^), uniformly on the contour C 
is uniformly 0(A^~"'^) on £2 and (3.20) implies 



^N^\0 is 0(A^ uniformly on the contour €2- Moreover, = 0(A^ so finally c\^\S,) 



-2s{z){m^^\z) - m(z)) = -{m^^\z) - m{z)f{z) + O ( sup |mj^^ - ) + 0{N-^). 



1^2 



Moreover, from the maximum principle for analytic functions, sup£2 |to^^ — to| < sup^^ jmj^'' — m| 
so the previous equation implies 



sup \mpf — m\ = O I sup \mpf 



We know that pi^^(s)ds converges weakly to p(s)ds (see [2]), so by (3.7) and (3.8) pi^'^^(s)ds 
converges weakly to /9(s)ds. On £1, z is at distance at least e from the support of both p'^'^^iyS^As 
and p(s)ds so, on £1, m)^ — m converges uniformly to 0. Together with the above equation, this 
implies that 

sup|m^^ - m| = O ^-^ 

By the maximum principle the same estimate holds outside Ci, in particular on £2, so equation 
(3.20) implies that for z inside £1 



-2siz){m^^\z) - m{z)) = -{m^^\z) - m{z)f + c^^\z) + O (^1) 



(3.21) 



Moreover, 



1 



1 



(z-A,)2 



< ^lmm^^(z) < ^|m^^(z) - m{z)\ + ^|Im to(z)| < -l-\m^^\z) - m(z)| + ^ (3.22) 



c 



for some constant c. We used the well-known fact that Im m is uniformly bounded on the upper 
half plane^. On the act A + 5 < E < B — 6 and |?7| < e. we have inf |s| > 0. Therefore (3.21) takes 
the form 



1 + 



(1. 



{m^^\z) - m{z)) = O ( \m''-\z) - m(z)|^ + -.k'^^^z) + 



From the hypothesis (3.15), if N~''-+'' < rj < s and A + S < E < B - 5, then 



|m^^ — m| < c|tojy 



(3.23) 



(3.24) 



/'DA 

for some c > and ejy as N 00. For large N, (3.24) implies that 



— m| > 1/c — 2£jv- Together with Im^-* — m\{E + ei) — >■ and the continuity of Itj^'^ — m| in 



m| < 2£jv or 



^This follows for example from properties of the Cauchy operator, see p 183 in [8]. 
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the upper half plane, this implies that Im)^ — < 2£jv and therefore Ini)^ — m| — > uniformly 
on N^^^'^ <ri<e, A + S<E<B — S. Consequently, using (3.23), this proves that there is a 
constant c > such that for any r] > N^^+"-, A + S < E < B - 6, 



\m^j^'\z) — m(z)\ < c 



1 



1 



(^) 



The same conclusion remains when substituting (resp. fc^'') by (resp. kj^) thanks to 

(3.7) and (3.8). 

□ 

To prove accuracy results for /i, the above Lemma 3.3 will be combined with the following one. 

Lemma 3.4 Let 5 < {B—A)/2 and E e [A+S, B~d] andO < rj < S/2. Define a function f ~ fE.rr' 
R — > M such that f{x) ~ 1 for x G (— oo,£' — 77], f{x) vanishes for x £ [E + 00), moreover 
|/'(a^)| < crj~^ and \ f"{x)\ < crj~^ , for some constant c. Let p be an arbitrary signed measure and 
let S{z) — J {z — x)~^p{x)dx be its Stieltjes transform. Assume that, for any x € \A-\-5/2,B — 5/2], 



\S{x + iy)\ < for ij < y < 1, and \ImS{x + iy)\ < for < y < rj. 
Assume moreover that J^p{X)dX = and that there is a real constant T such that 



(3.25) 



[-T.T]" 



|Ap(A)|dA < ^. 



Then for some constant C > 0, independent of N and E Cz [A + S, B — S], we have 

CU\logr,\ 



/B(A)p(A)dA 



< 



TV 



(3.26) 



Proof. Our starting point, relying on the Hclffcr-Sjostrand functional calculus, is formula 
(B.13) in [16]: 



fEiX)mdX 



<C 



yf'^{x)x{y)l^S{x + \y)dxAy (3.27) 



C {\fE{x)\ + \y\\f'E{x)\)\x'{y)\\S{x + iy)\dxAy, (3.28) 



for some universal C > 0, and where x is a smooth cutoff function with support in [—1, 1], with 

x{y) = 1 for \y\ < 1/2 and with bounded derivatives. 

Using (3.25) and (3.26), the support of x' being included in 1/2 < \y\ < 1, and the fact that 
is 0(77^^) on an interval of size 0(77), the term (3.28) is easily bounded by O (^). Concerning 

the right hand side of (3.27), following [16] we split it depending on < y < 77 and rj < y < I. 

Note that by symmetry we only need to consider positive y. The integral on the first integration 

regime is easily bounded by 



IO<y<ri 



2//ij(2;)x(y)Im5'(a: + iy)dxdy 



= O 



„ U \ 
yrj -——dxdy = O 



|2:-_E|<t;,0<i/<?) 



Ny 



'IT 
N 
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For the second integral, as and x are real, we can substitute Imm by m and use the analyticity 
of m when integrating by parts (first in x, then in y): 



< 



yfE{x)x{y)S{x + iy)dxdy 



ri<y 



yf'E{^)x{y)S'{x + iy)dxAy 



ri<y 



< 



dy{yx{y))f'E{x)S{x + iy)dxdy 
v<y 

'nf'E{x)x{ri)S{x + 'ni)dx 



This last integral is easily bounded by 0{U/N) using (3.25). Concerning the previous one, as 
f'E = 0(7,-1), „ ^1 < for non vanishing dy{yx{y)) - 0(1) and S{x + iy) = 0{U/(Ny)), 
this is bounded by 



O 



U f'dy 
Nj, y 



= O 



N 



concluding the proof. q 

As a corollary of Lemmas 3.3 and 3.4, we get the accuracy at scale 1/VN for the Afc's in the 
bulk. 

Corollary 3.5 For any a > and e > we have 

l7r-7.H0(iV-V2+.) 

uniformly in k ^ l_ceN, (1 — a)iV| where jlf^"^ and 7^ are defined by 



Pi = and J p{s)ds= — 



Proof. We will apply Lemma 3.4 to p = p — p[^'' with rj = N^^^"^^^ , and check the conditions 
on S = m — ruN ■ We denote z = x -\-\y. 

By the spectral gap inequality for the measure p, we get 



1 



\k=l 



Xk 



^^^^^ 



< 



(3.29) 



Together with Lemma 3.3, this implies that uniformly in N 1/2+^ < y < 1 and x £ [A + d/2, B — 
5/2], we have 



|mw(z)-m(z)| = 0(^^ + ^^ 
For < J/ < iV-i/2+'^, m is uniformly bounded and 



O 



Ny 



y ^ y lTamN{x + iy) 



7 L jP^^^Wdt 
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is an increasing function, so denoting yo ^ N ^/^+^ we have 

y |Im(TOAr(x+iy)-m(x+i2/))| < ?/o ImmAr(a;+iyo)+0(y) < yo\lm{mNix+iyo)-m{x+iyo))\+0{yQ). 
Therefore, for any < y < Af~^/^+^, 



|Im(mjv(a; + iy) ~ 'rn{x + iy))\ = O 



Ny 



Finally, the condition (3.26) with I] = □(TV^/^+e) ^nd with the choice of any T > max(|A|, \B\) + S 
follows from the large deviation estimate (3.5). 

Using the conclusion of Lemma 3.4 for functions ]e arid Je+ti defined in the same lemma, and 
subtracting both results, we get that uniformly in G [A + (5, i? — (5] , 



(pW(t)-p(t))dt 



0(iV 



-l/2+£ 



(3.30) 



so if 7^^^ e 8.B - S\, then |7^^^ - 7a;| = O {N-^I'^^^). This estimate holds uniformly 

in A; e [aiV, (1 — Q!)7V|: as a consequence of (3.30) and the smooth form (3.3) of p, for any 
k e |aiV, (1 — a)iV| and sufficiently large iV we have 7^^^^ € + (5, B ~ (5], for (5 > small enougii, 
concluding the proof. q 

Lemma 3.6 For any e > there exists ci,C2,e' positive constants such that for any iV^/s+e <^ 
j <N - 7V3/5+'^^ we have 

Proof. We will assume that j < N/2 in the following, i.e. we will estimate the accuracy near 
the edge A, the proof close to the other edge B being analogous. We follow the notations used in 
Corollary 3.5 and Lemma 3.4. For e [A ~ S, A + 6] Li [B - S, B + 6], we have inf \f\{z) > 
Therefore we can divide — 2/(z) on both side of (3.21) to have 



1 + 



/ 1 



ViVr/3/2 



(™jv (2) " "1(2:)) = — 



By (3.29), (3.7) and (3.8) we can bound the variance term by 



(3.31) 



(3.32) 



for r] > N ^'^ for example. Following the same continuity argument in the proof of Lemma 3.3, we 
obtain that for any e > 



\lm{mpf — m){x + i?7)| = O 



provided that 



1 



7V?y3/2 



< 



< 1. 



(3.33) 
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We can now follow the argument in the proof of Corollary 3.5 so that (3.25) holds with U — iV'^/^. 
Since the condition (3.26) is easy to verify, we thus have 

' Mt)[pi^Ht) - p{t)]dt <^^, v-N-y\ 

where /e is defined in Lemma 3.4. This proves that 



In particular, 



N 7V2/5 J.ao ^'/^ 

and we have, by definition of ji, that 

Similarly, the reverse inequality holds and we have 

7j_Ar3/5 + . - 277 < 7]"'' < 7j + 7V3/5 + e + 277. 

Since p{t)dt ^ {E - Af/^, for j > N^/^+^ we have 



|7j-Ar3/5+» - 7i 



Moreover, by (3.12), Xj is concentrated around Ep(Aj) at scale N ^Z^, so |E^(Aj) — 7j^''| = 
0(A^~i/2), concluding the proof of the lemma. q 

3.2 The locally constrained measures 

In this section some arbitrary e, a > are fixed. Let be a continuous nonnegative function with 
6* = on [-1, 1] and 6*" > 1 for |a;| > 1. We can take for example 0{x) = {x-l)'^lx>i + {x+l)H^^^i 
in the following. 

Definition 3.7 For a given k S faN, (1 — q;)A'^] and any integer 1 < M < aN , we denote 
j(fc,A/) = |fc _ ^_ and iM = \l(k,M) \ ^ 2M + 1. Moreover, let 



We define the probability measure 

do.C^'^^) ^e-^'"-"" dp, (3.34) 

where Z = Z^^i^k.M) . The measure uj^^'^^^ will be referred to as locally constrained transform of p, 
around k, with width M . The dependence of the measure on e will be suppressed in the notation. 
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We will also frequently use the following notation for block averages in any sequence xi,X2, ■ ■ ■ 

4''':=^ E (3-35) 

The reason for introducing these locally constrained measures is that they improve the convexity 
in J^'^^^^) up to a common shift, as explained in the following lemma. 

Lemma 3.8 Consider the previously defined probability measure 

Z Z 

where we denote 

^1 := ^^^^^^^ - f E log|A.-A,|, 

'^^■■=-^ E log|A.-A,| + |E^(^«) 

where jC^^^^) is the set of pairs of points i < j in |1,-/V| such that i or j is not in and 
Hi = HiiXk-M, • ■ • , Afc+7\/). Then V^'H2 > and denoting v = {vi)^^i{k.M) , we also have 

v*{^^ni)v>!^— E i^—^j)'- (3-36) 

Proof. Since V is convex, to prove the convexity of , it suffices to prove it for the Coulomb 
interaction terms; this relies on the calculation, for any u £ M.^ , 

Moreover, for any v S R**' , a similar calculation yields 

MVHi).>^ E jx—x;)2+^^ E (-^--.)^(^ — 

(3.37) 

From our definition of 0, 

1 iV2-2%„ /A^i-^A, - A,)\ ^ iV2~2e 

> — . 



{Xi — Xj)^ i\j \ IM / I'M 

which implies 

jYl-2e 

v*{'^^ni)v>p—^ E i.'"^~^Jf■ (3-38) 



i<j,i, 



which completes the proof of (3.36), noting that the above factor 1/2 comes from the strict ordering 
of i and j indexes in (3.38). q 
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The above convexity bound, associated with the foUowing local criterion for the logarithmic 
Sobolev inequality, will yield a strong concentration for "^i^i under o;'^'^'*^-', if = 0. 

This lemma is a consequence of the Brascamp-Lieb inequality [7] . Notice that the original inequality 
applied only to measures on K^, but a mollifying argument in Lemma 4.4 of [17] has extended it 
to the measures on the simplex {Ai < A2 < . . . < \n} considered in this paper. 

Lemma 3.9 Decompose the coordinates A = (Ai, . . . , A^r) of a point in = R™ x M^"™ as \ ^ 
(x, y), where x € K™, y G R^"". Let oj = ^e"^'" be a probability measure on R^ = R" x R^"™ 
such that % ~ Hi + ^^2; with Hi ~ 'Hi(x) depending only on the x variables and 'H2 = '^2(2;, y) 
depending on all coordinates. Assume that, for any A G R^, V^'H2(A) > 0. Assume moreover that 
'Hi(x) is independent of xi + . . . + Xm, i-c. X^i'^i di'Hi{x) = and that for any x,v G R™, 

^ m 

V*{\7^ni{x))v > — ^\v^- Vj\^ (3.39) 

with some positive ^ > 0. Then for any function of the form /(A) = "ViXi), where Wi = 

and F : W W is any smooth function, we have 

I /^log/^do. - (^J fduj^ log f^du^ ^ ^ / l^/I'd^- (3-40) 
Proof. In the space R™ we introduce new coordinates {z,w) = M*{xi, . . . ,Xm) with z = 

(zi,...,z,„_i) eR™-i , 

w:=m-i/2^Xi, (3.41) 

i 

and M is an orthogonal matrix. Since 'Hi{x) is independent of I'l + . . . + Xm, we can define 
'Hi{z) := Hiix). Similarly, the function /(A) = Fi^^^ ViXi) with J2i = depends only on the 
z coordinates, i.e. it can be written as g{z) = /(A). Hence we can rewrite 

fHogfdu;= [ logged,,, [ fdiu^f g^dv, 

where dv = v{z)dz with 



Z 



xR"- 



Introduce the variable q = (w, y) G R x R ™ and denote by 'Hqq,'Hzq,'Hzz the matrices of second 
partial derivatives. As "^2 is convex, the Brascamp-Lieb inequality yields 



nzz>^ [ e-^«(^'^) IHzz - Hzq[Hqq\-^n 



dwdy. 



Since Hi is independent of q, we have, by assumption of the positivity of the Hessian of ■H2, that 
for any q fixed, 

(H2).. - Uzq\nqqY^Uzq = (^2).. " (^2 ) .9 [(^2 )99] " ^ (^2 ) > 0. (3.42) 

Thus we have, for any u G R™^"'^, that 

u^UzzU > u*{Hi)zzU = u*M*{ni),,Mu > ^ y [(Mu), - (Mw)jf , (3.43) 
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where M denotes the first m — 1 columns of M . Since the last column of M is parallel with 
(1,1,...,!) € R™ and M is an orthogonal matrix, we have X^il-^^"")* ^ ^ ^^^^ 

^ m m m — 1 

5] [(Mh), - {Mu),f = 2e5][(M^.),]' - 2e 5] (3.44) 

2,_7 — 1 z— 1 i—1 

Hence the measure v ~ exp(— A^H) is log-concave with a lower bound 2NS^ on the Hessian of NT-L^ 
and we can apply the Bakry-Emery argument to prove the logarithmic Sobolev inequality for v. 
Without loss of generality we can assume that j pduj = J g^di' = 1. Therefore, we have 

/ fHogfdu;^ [ g^\ogg^Av<^ [ |V,5pdi. / |V,/pdc., (3.45) 

where we have used the orthogonality of M to show that jV^^p = IV^/P- This proves the estimate 
(3.40). □ 

It is now immediate, from Lemma 3.8, Lemma 3.9 and Hcrbst's lemma, that the following 
concentration holds. 

Corollary 3.10 For any function f{{\i, i G /f*^^^^)}) = 'J2ii''.'^n ''^i^i with Vi — Q we have 

\ 4iM|wr 

Choosing Vj ~ ^Vj+i = 1 and all other u^'s being zero, this corollary shows that the particle 
differences \j — concentrate around their mean with respect to the o;''^'^^-' measure. By 

choosing e small and M almost order one, we obtain concentration almost up to the optimal scale 
If we can justify that the measures oj^'^'^^-' and fi are very close (in a sense to be defined), we 
will have concentration of differences at the optimal scale for fi. We will then separately show, by 
using the loop equation, that accuracy will hold at the same scale as well. This is the purpose of 
the next subsection, through an inductive argument. 



3.3 The induction 

The purpose of this paragraph is to prove the following proposition: if accuracy holds at scale 
A/'-i+°, it holds also at scale TV-^+t". 

Proposition 3.11 Assume that for some a G (0, 1) the following property holds: for any a, e > 0, 
there are constants 5, ci, C2 > such that for any N > 1 and k £ {ctN, (1 — Q;)Af], 

P,. (|Afc - 7fc| > N-'+''+') < cie-^^^'. (3.46) 

Then the same property holds also replacing a by 3a/4; for any a,e > 0, there are constants 
(5, Ci, C2 > such that for any N > 1 and k £ \aN, (1 — a)A^], we have 

P,, (|Afc - 7fc| > N-^+i^+-) < cie-^^^'. 
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Proof of Theorem 3.1. This is an immediate consequence of the initial estimate, Proposition 
3.2, and iterations of Proposition 3.11. q 
Two steps are required in the proof of the above Proposition 3.11. First we will prove that 
concentration holds at the smaller scale Af~^+t. 

Proposition 3.12 Assume that (3.4S) holds. Then for any a > and s > 0, there are constants 
Ci,C2,6 > such that for any N > 1 and k G faN, (1 — Q;)Af], 



|Afe-E^(A,.)| > j <cie 

The above step builds on the locally constrained measures of the previous subsection. Then, 
knowing this better concentration, the accuracy can be improved to the scale N~^^^ . 

Proposition 3.13 Assume that (3.4-6) holds. Then for any a > and e > 0, there is a constant 
c > such that for any N >1 and k e faN, (1 - a)N}, 



{N) 

Ik - 7fe 



< c- 



+e 



N 



Proposition 3.11 is an immediate consequence of the last two propositions. The proofs of these 
two propositions are postponed to the end of this section, after the following necessary series of 
lemmas. 

Lemma 3.14 Take any e > and a > 0. There are constants ci,C2 > such that for any 
N >\, any integers 1 < Mi <M< aN , any k e {aN, (1 - a)Nj, and wC"'-^^) from Definition 3.7 
associated with fc, M, e, 



Proof. Note that a[!*^^' — x\^'^^ is of type J2ii''.'^n ''^i^^i with = ^''^^ 

1 ^ ^ Mi + l 1^1 / + 1 

Together with Corollary 3.10, this concludes the proof. q 

Lemma 3.15 Assume that for /i accuracy holds at scale N'^^"", i.e. (3.46). Take arbitrary 
a,e > 0. There exist constants ci,C2,S > such that for any N > 1, for any integer M satisfying 
< M < aN/2, for any k G faiV, (1 - a)N\ and for any j G |1, Nj, we have 

|E,,(A,) -E^(..«)(A,)| < cie-=^^', 

where the measure o;'^'^'"'^^-' is defined in (3.34) from Definition 3.7 with parameters k,M,e. 

Proof. First, the total variation norm is bounded by the square root of the entropy, and by (3.8) 
the particles are bounded with very high probability, so we have 




|E^(A,) - E^(.,M, (A,)| < C'v/^(M|a.(Mf)) + o(e-^^) 
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for some c, C > independent of k,j. For the measures we are interested in, using the logarithmic 
Sobolev inequaUty for /i, we have for some c, C > 

Now, as d"{x) = if |x| < 1 and 9'{x)'^ < 4a;^, for some new and universal constants c, C > 



< CN" 



IE^((A 



I Afc + M —^k-AI |> 

1/2 



k+M 



Afc-A/)^) ( |Afe+7i/ — Afc_j\/| > 



1/2 



(3.47) 



This moment of order 4 is polynomially bounded, for example just by concentration of order N ^1"^ 
for all Aj's under /i. Concerning the above probability, as |7fc+A/— 7fe-M| = 0(il//iV), for sufficiently 
large N if |Afc+A/ - Afc_A./| > -^¥r7 then either |Afe+M - 7fe+M| > M/N'^'^ or |Afc_A/ - 7fc_A/| > 
M /N^~^ . But accuracy holds at scale iV~^+° < M/N, so both previous events have exponentially 
small probabilities, uniformly in k. Indeed, one has k — M, k + M E laN/2, (1 — a/2)N^ and by 
(3.46)there are constants S, ci, C2 > such that for any > 1 and k £ laN/2, (1 — a/2)N^, 



{\Xk-jk\>M/N^~') <cie 



This concludes the proof. 



□ 



Lemma 3.16 Assume that for ^ accuracy and concentration hold at scale N ^+<^ , Take arbitrary 
a,e > 0. There are constants ci, C2, 5 > such that for any N > 1, any integers N"" < M < aN , 
1 < Ml < M, and k e paN, (1 - 2a)Nl, 



x[^- - Ai: 



[A/] 



[Ml 



Afl) 



> 



"/V V 



< cie 



Proof. Consider the measure uj^'^'^'^') associated with the choice e = e/2. First note that, by 
Lemma 3.15, 



E 



A 



[A/] 



E, 



[Ml 



< ce 



for some coefficients c,C',di, uniformly in N, k, M, Mi. As a consequence of this exponentially 
small difference of expectations, the probability bound to prove is equivalent to the existence of 
ci , C2 , i5 > such that 



(A) < cie- 



A 



[M{ 



A 



[A/] 



E 



> 



[M 



N V Ml 



with the same uniformity requirements. By Lemma 3.14, we know that there are such constants 
with 

P^(..M, (A) <cie-=^^', 

so the proof will be complete if we can prove that |P(^(fc,j\/) {A) — P^j. [A) \ is uniformly exponentially 
small. By the total variation/entropy inequality we have: 



|P^(.,M, {A) - P^ (A) I < / Idc.^'^'^^) - d/i| < ^2S{i^^\J*^). 



(3.48) 
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This entropy was shown to be exponentially small in the proof of Lemma 3.15, see equation (3.47). 
□ 

Lemma 3.17 Assume that for ji accuracy and concentration hold at scale N^^^"" . For any e > 
and a > 0, there are constants ci, C2,S > such that for any N > I and k e l2aN, (1 — 2a)N'l, 



> 



N 



< cie 



Note that in this lemma and its proof, for non-integer AI we still write AJ, for AJ.*^ ^ , where [MJ 
means the lower integer part of M . 
Proof. Note first that 



A,~Ar^l-E,(A,-Ar') < A,-Ari-E,(A.-Ar') 



Ar-Ari-E,(Ar-Ari) 



By the choice Mi — 1, AI — iV" in Lemma 3.16, the probability that the first term is greater than 
— is exponentially small, uniformly in fc, as desired. Concerning the second term, given some 
r > and q N defined by 1 — r < a + < 1, it is bounded by 



9-1 

i=0 



k ~ \ 



V ~ -^k 



[a+qr] . [aN] 



By Lemma 3.16, for any e > 0, each one of these 9 + 1 terms has an exponentially small probability 
of being greater than ^ ^ . Consequently, choosing any e and r (and therefore q) such that 

□ 



e + ^ < a/2 concludes the proof. 



Proof of Proposition 3.12. We just need to write 



|Afc-E^(Afc)| < 



Afc-Ai"~l-E,(A,-Ari) 



lE.(Ar') 



By Lemma 3.17, the first term has exponentially small probability to be greater than " ^' . By 
the logarithmic Sobolev inequality for /i with constant 0(1) (see (2.2)), the second term has an 
even better concentration, at scale 1/A^: 



(iaL" 



-E.(Ari)l>|)<Ce- 



This concludes the proof of the proposition. 



□ 



Proof of Proposition 3.13. Thanks to Lemmas 3.3, 3.4, and reproducing the proof of Corollary 
3.5, we know it is sufficient to prove that for any 5 > 



1 

iV2 



var 
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goes uniformly to where z = E + irj, EE [A + S, B — 6] and 77 > iV^^+T^+'^. 

Let io be the index in [0, A^] such that the typical position 7,^ is the closest to E. Define the 
indexes of particles close to E, far from E and in the edge as 

Int = {i: \i~io\ < iV"+^}, 

Ext = {i: \i-io\ > iV"+% i G [aiV, (1 - a)A^l}, 
Edg = {t:i^ laN, (1 - a)Nj}, 

where a is small enough such that 



6 S 
jaN<A+-<A + S<E<B-6<B~-< 7(i_„)jv- 



(3.49) 



We choose ak = E^(Afe) in the following equations. Then 



1 

iV2 



z - Ai 



fcGEdg 



-A, 



fee Ext 

The edge term is bounded by 



z - Xk z - ak 



z- ak 
C 

N2 



fcGint 



z - Xk z - ak 



(3.50) 



— > 

Edg 



1 



1 



z - Xk z- ak 



Edg 



E-Xk\<- 



c 

NS- 



J2E,{\Xk^akn 



Edg 



(3.51) 

From the condition (3.49) and the large deviation estimate (3.5), the above probability is expo- 
nentially small. Moreover, the above moments are 0(1 /N) by the spectral gap inequality for 
/i, see e.g. equation (3.12). Hence the edge term goes to uniformly. 

Using the accuracy at scale N~^~^'^ and the concentration at scale Af~i+°/2 (Proposition 3.12), 
the second term in (3.50) is bounded, up to constants, for some ci, C2, (5 > by 



Ar2 



E 



Xia+k — Oii„ + k 



(±\2 



Cie 



E 

, fc>Af°+' 



\Xto+k - "io+fcl^ 



E 



cie 



'C2N'' 



In particular, it converges uniformly to 0. For the third term, for some c > it is less than 
1 



■E 



\Int I j ' \lnt / 



< C 



{N' 



a+6\2 



2e 



This goes to if 77 1+ ^4" + 2 , concluding the proof. 



□ 
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4 Local equilibrium measure 
4.1 Construction of the local measure 

Let < K < 1/2. Choose g G [k, 1 — k] and set L = [Nq] (integer part). Fix an integer K with 
K < {N — L)/2, in fact we wih always assume that K depends on. N a.s K = TV*^ with fc < 1. We 
will study the local spacing statistics of K consecutive particles 

{A, : J el}, I = Il ■■= IL + 1,L + Kl 

These particles are typically located near Eg determined by the relation 




p{t)dt = q. 



Note that |7l - Eq\ < C/N. 

We will distinguish the inside and outside particles by renaming them as 

(Ai, A2, . . .,\n) ■■= (2/1, ■ ..yL,XL+i, ■ ..XL+K,yL+K+i, ■■■Vn) e S^^^ (4.1) 

but note that they keep their original indices. The notation S^^^^ refers to the simplex {z : zi < 
2:2 < . . . < zjv} in ■ In short we will write 

■>c= {xl+i,---Xl+k), and y = {vi, ■ ■ ■yL,yL+K+i, ■ ■ -Vn), 

all in increasing order, i.e. x G S'^' and y S ^(n-k)^ .^jjj^ refer to the y's as the external 
points and to the x's as internal points. 

We will fix the external points (often called as boundary conditions) and study the conditional 
measures on the internal points. We consider the parameters L and K fixed and we will not 
indicate them in the notation. We first define the local equilibrium measure on x with boundary 
condition y by 

(4.2) 

where u is the density of /iy. Note that for any fixed y S xj lies in the interval 

[yL,yL+K+i]- 

Given the classical locations, 7 = {71, 72, . . . , 77v} with respect to the /^-measure, we define the 
relaxation measure fi^'^ = fi'^ by 

dM^ := ^e^^Q^dA^, Q-(x)=^QJ(.T,), QJ(:,) = -l(a: - 7^)2, (4.3) 
Z„T ^ — ' ■' ■'It 

where is chosen such that /i is a probability measure. Here < r < 1 is a parameter which may 
even depend on y, i.e., r = r(y) is allowed. Note that an artificial quadratic confinement has been 
added to the equilibrium measure. We define the local relaxation measure /iy to be conditional 
measure of fj,'^ . 

Define the Dyson Brownian motion reversible with respect to p,y , by the Dirichlet form 

^M;(/) = E^/(^'/)'dM;, (4.4) 
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Aiy(dx) = Uy{x)dx, '"y(x) := u{y,x) 



J u(y,x)dx 



where di = dx^- The Hamiltonian "Hy of the measure ^y(dx) ~ exp(— iVHy) is given by 

ni{^) = Y,^Vy{x,)-^Y. ^og\x,-x,\+Y,Ql{x,), (4.5) 

yy(x) = T/(x)-l^log|x-y,|. (4.6) 

We now define the set of good boundary configurations with a parameter eg > and a parameter 
(5 = 5(A^) > that in the apphcations may depend on A^: 

e^.so = ^ := {y e S^^""^^ : \yk - lk\ < 5. "ike IN K/2,LlyjlL + K + l, Nil -n/2)l (4.7) 

and \yk-lk\ < 1, Vfce |1,A^1, 

and E^^ {xj - 7j )2 < for all j G [L + 1, L + iiT] 

and yL - yL-i > cx.p{-N^"), yL+K+2 - Vl+k+i > exp(-Af^'')|. 

First we show that the good configurations have overwhelmingly large probability 

Lemma 4.1 For any £o > and for any choice 5 = N^"^ with d G (1 — fc, 1), there is an e' > 
depending on d such that 

Pm(S') < Ce-"^'' + Ce-"^'" . (4.8) 

Proof. We have proved in Theorem 3.1 that for any choice S = TV"'' with d € (0,1) the 
probability that the first condition in (4.7) is violated is bounded by Cexp(— cA'^^ ) with some 
e' > depending on d. Similarly, the second condition is violated with an analogous very small 
probability by (3.5). To check the probability to violate the third requirement in the definition of 
Q, we use that 

P^{e^^(x, - 7,)2 > 5'} < P,,{P^,{|x, - 7,1 > 5/2} > SS^a} + Cexp(-ciV-') 

< C6-^E^F^,^{\xj - 7jl > 5/2} +Ccxp{-cN'') 

< C5-^¥^{\x,~^,\ > 5/2} < cie-^^^^', (4.9) 

since for y satisfying the first two conditions of (4.7) we have 

^f^A^J - Ij)' < -574 + - 7,1 > S/2} 

as Xj — jj < y^+if+i — 7j < (5 + "fL+K+i ^ 7i < 1 and also a similar lower bound holds. 
Finally, we show that 



{yL-yL-i<cxp{-N'«)) < Ce 



and a similar bound holds for the other condition in the fourth line of (4.7). For simplicity of the 
presentation and to avoid introducing new notations, we will actually prove 



.(j/L+i-J/L <cxp(-iV^")) <Ce 
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from which the previous inequality follows just by shifting the indices. With the events 
A := {yi+i - 2/i < exp(-iV^'')} , Q := {vl+k+i ~ Vl < 2a}, 

we write 

P^yl) =E4i(17)P^^(A)] +P^(r!^). (4.10) 

Choosing a ~ CqK/N with a sufficiently large fixed constant Co Theorem 3.1 and 6 <^ K/N 
guarantee that P^i(il^) is subcxponentially small. 
We will prove that 

[xL+i -VL< N~\) < Cvr (4.11) 
for any r G (0, 1). The constant depends on V , more precisely 

Cv = C + Csu^{\V'{x)\ : x^[yL,yL+K+i]]. (4.12) 

From (4.11) the necessary subexponential estimate on the first term in (4.10) follows by choosing 
r = A^-2gxp(-7V'^o). 

To prove (4.11), on the set Q. we can shift the measure such that = —yh+K+i and denote 
a :~ — J/L- Then we have 



J J —a-\-aip .j 



where we set Wj := (1 — ip)~^XL+j, dx = dx^+i . . . d-T^+x and dw = dwi . . . dwK- By definition, 

i<L i>L+K+l 

i<L i>L+K+l 

Note that we only used that V is a C^-function with bounded derivative in performing a Taylor 
expansion and using that Wj < a is finite. Hence 



Z 

with 



1 / . . . r dx n (2;. - X,fe-^^ > (1 - ^)NK+CK-^-CvNK^ 

^ J J —a-\-aip ijGi 

i<3 



z := r dw n ( 



Therefore the /iy-probability of i/l+i yL ~ — yL > a{\ — ip) can be estimated by 

P^A^L+i >-a + ^a) > (1 - ^)NK+CK-^-c.NK^ >l-{Cv + C)NKip 
by using K < N. Choosing (p = N^'^r/a and recalling that a ~ K/N, we arrive at (4.11). q 
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Proposition 4.2 Let > be fixed. For any smooth, compactly supported function G 
we have 



lim 



N- 



0, 



provided 



1 



!i '''(y) 2f /or any y ^ Q 



holds for the function r = T(y) mi/i some constant t = tn such that 



(4.13) 
(4.14) 

(4.15) 



We remark that, with a shght abuse of notation, tlie last term, i = L + K in tlic sum involving 
the non-existing Xi+i = xl+k+i is defined to be zero. We also point out that the notation E^E^^ 
means that the law of y is given by ^ in the first expectation and we are using the measure fiy in 



the second one. Of course, we have E^ = E^E^^. 



Proof. For any configuration y, any t (may depend on y) and for any smooth function G with 
compact support, we have 



[E^^-E„j]-i^G(iV(x,-x,+i 

iei 

Here we also introduced the notations 
D{fi I w 



< C 



K 



D{^ly\,ll))"\ce-^''"' ^ S{iiy\yil), (4.16) 



1 

27V 



2N 



(4.17) 



and 



S{^Ji\u^) := j log(^)dM 



for any probability measures The estimate (4.16) follows from our the local relaxation to 

equilibrium argument that in this form first appeared in Theorem 4.3 of [15]. We will neglect the 
exponentially small entropy term since it can be estimated by the Dirichlet form, i.e. by the first 
term as long as r > N~'-^ . 
We thus obtain 



1 / N /NV/"^ \l/2 

E4E^^-E^j]-^g(7V(x,-x,+i)) <c(-^E4ler(y)Z?(Aiy|A4)]) +Ce 

iei 



-cN" 



(4.18) 



To obtain the estimate (4.18) we separated good and bad configurations; we used (4.16) for y G G- 
On the complement Q'^ we just used the trivial estimate on G, and this yields the subexponentially 
small second term. 

Assuming (4.14), we have 



K 



< — < N-^ (4.19) 



by (4.15). Inserting this estimate into (4.18) we completed the proof of the proposition. 



□ 
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4.2 Matching the boundary conditions 

Suppose wc have measures a and /i with potentials W and V given by (2.1) with densities p — pv 
and pwi respectively. For our purpose W{x) ~ , i.e, a is the Gaussian /3-ensemble and pw{i) = 
^ \/[4 — ig the Wigner semicircle law. Let the sequence 7j be the classical location for /i and 
the sequence dj be the classical locations for u. 

We will match the boundary conditions for the local measure on Jy := [j/l, around 
Eq with those of the a measure. For definiteness we choose the interval J' = [O^i ,9li+k+i] with 
L' = i(iV — A' — 1) as our reference interval. Note that J' is symmetric to the origin. The local 
density pv{Eq) at the point Eq we look at may be different from the density pw[^) at the origin. 
Thus the typical length of Jy, which is ^l+k+i — 7l ~ [pv{Eq)]~^N~^ , may not be close to the 
length of J' which is very close to [pvy(0)]~^A^^^ = ttN"^, so we will have to rescale the a measure 
by a factor 

_ pyjEq) 

' ~ pw{0) ■ 

In fact, we need to match not only the interval of classical locations 7 with J', but the exact 
interval ly. This requires a y-dependendent scaling factor s = s(y). 

From now on we assume that y is a good boundary condition with a parameter 5 that satisfies 

SN 

— ^ 0. (4.20) 

We can shift the coordinates so that 

-yL = yL+K+i. (4.21) 

Since our observable is translationally invariant, we will not track the translation and we assume 
that (4.21) holds. We define 

.(y) := ^ = sq ^. (4.22) 



VL VL+K+l IL 



We have 



\s{y) 



VL IL 



SN 

< C-— (4.23) 
K 



Ql' « -\ewm-'^. IL « -[Qv{Eq)r'^, VL « AQv{Eq)]-'^, (4.24) 

by using y G ^ and (4.20). Similar formulas hold for O^ij^k+i^ 1l+k+i and hl+k+i at the upper 
edge of the interval. Here the A k B is understood in the sense that the approximation error at 
most of order (K/N)"^, recalling that K = o{N). 

For simplicity of the presentation, we can first shift the original /i-ensemble such that Eq = 0. 
Second, we can perform an initial rescaling of the Gaussian /3-ensemble so that Sq = 1. 

Lemma 4.3 Assuming Eq ~ 0, Sq ~ 1, we have 

llL+j - 0L'+j\ <cj^ + CS, \j\ < ^N^. (4.25) 
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Proof. The classical locations jL+j and Ol'+j are given by the equation 

Pv{^) = j^, Pwi^) = j^. (4.26) 

We will use the approximations 

Pv{x)=pvi0) + O{x), pw{x)=pwiO) + 0{x) (4.27) 
for small x (to stay away from the spectral edge). Since \yj ~ "fj\ < i5, we have 

-^=J pv{x)dx^ j pv{x)Ax + 0{5)^ pv{Q){vL+K+i-VL) + Oij^\^-0{5). 

'"'^ (4.28) 

Similarly, 

= Pw{x)Ax = pw{mL'+K+i - 6v) + O j . (4.29) 

Since —yh = Vl+k+i = ^^L'/s = Qu+k+iI s which is comparable with K/N by (4.24), and since 
|s - 1| < from (4.23), we have 

|pw(0)-py(0)|<^ + ^. (4.30) 
From (4.26), (4.27) and (4.24) we get 

■2 ■ -2 

which combining with (4.30) and pwi^) ^ c gives 

■2 

- 7L = Ol'+, - e^, + O(^) + 0(5). 

Since 7^ ~ ^l', this completes the proof of the lemma. q 



4.3 Rescaling of the reference problem 

Throughout this section we fix a good boundary configuration, y € G and a number r(y) depending 
on this configuration and satisfying (4.14). We will approximate the local relaxation measure Py 
on [yL, Vl+k+i] by a fixed reference measure. 

Given the collection of classical locations 9j corresponding to the Gaussian potential W{x) = x^ 
we define a reference local relaxation measure ag via the Hamiltonian 



Til 



«'W = E[f^'-|^El°8>''-^^l]-| E log|x,-x.| + ^^(x.-^.)', (4.31) 

iel' j^/' i.jei' is/' 

i<i 

on the set [9l' ,Ol'+k+i\ where /' := |L' + 1,L' + K\. Note that if a is the equilibrium measure 
given by (2.1) corresponding to W and cr"^ denotes the corresponding relaxation measure 

da' := -^e-^Q\la, g^(x) = ^ Q;(x,), gj(x) = ^{x - O.f, (4.32) 
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defined analogously to (4.3), then is the conditional measure of cr"^ under the condition that the 
outside points are exactly at their classical locations, i.e. Xj = 9j, j ^ I'. 

We make three simplifications in the presentation. First, as already in Section 4.2, we assume 
that both the configuration space [yL^yL+K+i] for the original measure /Xy and the configuration 
space ^L'+if+i] of the reference measure <7g are symmetric around the origin; this can be 
achieved by an irrelevant shift. Second, we assumed Sq = 1, which can be achieved by an irrelevant 
rescaling of W . Finally, we will set L' — L. This last assumption expresses an irrelevant shift in 
the labelling of one of the ensembles. Strictly speaking, shifting would mean that the original set 
of particles indices |l,iV| gets shifted as well. However, in our argument this shift does not play 
any active role; the only information we use about the set of indices is that L is macroscopically 
separated from its boundary and that its cardinality is A^. 

We now rescale the measure crj from [9l,0l+k+i] = to [yL,yL+K+i] = [yL^-yh] by 

the factor s = s(y) defined in (4.22) (note that y^, < 0). With the rcscalcd boundary conditions 
9j — !■ O'j := 9j / s, we define the reference local relaxation measure, or reference measure in short, to 
be 

^^-^''^■^^^^dx, (4.33) 
a measure on the set [yL,yL+K+i\ with Hamiltonian 



n 



"'^(x) = E - ^ E log - ^./^i] - 1^ E i°g 1^^- - + E(-^' - ^»/^)'- (4.34) 



The rescalcd potential associated with this Hamiltonian is Ws{x) = s^x^. 
For any smooth function G with compact support, we have 

E4'4E^(^(^^-^-'+i)) =^ / ''^'dxe-^«^-°('') }^J2G[Nix,-x,+,)), (4.35) 

where dx stands for the ii'-dimensional integral Jj^ diCi+i . . . dx^+x and Z"^'^'" is the nor- 
malization factor. Let Xj — Wj/s, then the right side becomes 



dw e-^-^?(-)-i Eg( ^("' ^ "'+^^ ) = Ki^T^g C^"' ; ""''^ ) (4.36) 



" iG/ ■ " " iG/ 

= E 



where we renamed the w-variables to x-variables in the first step and in the second step we have 
used that 

g[n{x, - a;,+i)/s) - g[n{x, ~ x-,+i)) \<C\1- s\\\G'\\oo 

by Taylor expansion and from the fact that G is compactly supported. Clearly, the difference 
vanishes as long as s — > 1. Thus we are free to scale the measure with factor converging to 1. The 
condition s ^ 1 will be guaranteed by (4.23). 
Our main result is the following theorem. 
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Theorem 4.4 LetO < ip < ^. Fix K ^ N'' , S ^ N-'^, f = TV"* with d=l-ip,t^ 2d-l-ip^ 
1 — 3(^ and k ~ '^f, "in particular such that (4.15), (4.20) are satisfied. Then 



zG/ is/ 



(4.37) 



as — > oo /or any smooth and compactly supported test function G. Here the law ofy is given by 
fi in the expectation. 



Proof. From the rescahng estimates, (4.35)-(4.36), it suffices to prove that 
E^[E^j -E^...(,)]-i^G(A^(a;,-a^,+i)) ^0 



(4.38) 



as N ^ oo. Notice that after the rescahng both measures fXy and cr^'^^'''' hve on the same interval 
[vl, Vl+k+i]- In Proposition 4.2 we already showed that 



E, 



, [E^r - E^,/.,,>2] 1 ^ g(7V(x, - x,;+i)) ^ 0, 



(4.39) 



since (4.23) with Sg = 1 and 6N/K — !• guarantee that r(y) := f/s(y)^ satisfies (4.14). Thus the 
limit (4.38) will follow from the following Proposition that we will prove in Sections 5: 

Proposition 4.5 Under the assumptions of Theorem 4-4j we have 



EJE 



E 



This completes the proof of Theorem 4.4. 



(4.40) 



□ 



Proof of Theorem 2.1. Finally, combining Theorem 4.4 with Proposition 4.2 and noticing 
that (4.15) is satisfied since t = 2(i — 1 — we have 



E^-J, ^ G(iV(a:, - x,+i)) - E,. ^ G(iV(x, - x.+i) 



iei 



-+ 



(4.41) 



as — > oo. This holds for K = A^*^ with any < fc < | by selecting a suitable ip in Theorem 4.4. 
However, the measure ag is independent of V, the only information we used was that the local 
density matches. So we obtain that any two measures fJ-p^v and fip^w have the same local gap 
statistics assuming that the local densities of the two ensembles coincide. q 



5 Comparison with the reference problem 

In this section we prove Proposition 4.5. On the set y g CJ'^ with subexponentially small probabality 
(4.8) a trivial estimate on G suffices. For the sequel we therefore assume that y G G and we set 
T(y) := f/s(y)^ which clearly satisfies (4.14). In the first step we will soften the boundary condition 
y for the measure local relaxation measure /Xy. 
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5.1 Regularizing the boundary conditions 

We know that the boundary condition y G is regularly spaced on the scale S — N~'^, but it 
does not exclude that NS — N^^"^ ^ 1 points of the coUetion y pile up near the edges of the 
interval [y^, ijl^k+i]- This would substantially influence the local relaxation measure /iy near the 
corresponding edge inside [j//,, j/^+a'+i]- We therefore first replace the boundary conditions near 
the edges by the regularly spaced ones given by 9' = 0/s{y). This change will be controlled only 
in the entropy sense. The local relaxation measure with regularized boundary conditions will then 
be compared with the reference measure in the stronger Dirichlet form sense. 
Set a parameter 

B ^ N'' with l + ip-d<b<k, (5.1) 

in particular SN <^ B <^ K. Given a boundary condition y G ^, we define a new boundary 
condition ^ {vf ■ i ^ 1} as 

{max{6'^, yL-is} for L-AB <i<L 
y, i<L-4B, or i>L + K + AB (5.2) 

mm{ei, yL+K+4B} for L + K + l<i<L + K + 4B, 

i.e., we replace at most AB boundary conditions yt with the rescaled classical ones 9[ = Oi/s{y) 
near the edges of the interval [yL,yL+K+i\ = [6*^, ^i+^+i]- Note that the configuration space is 
unchanged. We have 

VL-AB < iL-iB + 5 < 9'^_^s + CB^N-^ + CS< 9'^_^s, 

where we used that y £ ^ in the first step and (4.25) in the second. In the last inequality we used 
that 0'j^_2B ^ ^'l-ab — cBN^^ (by regular spacing) and the definition of B from (5.1). Thus we 
obtain 

yf^el, L-2B<t<L, (5.3) 

and similarly at the upper edge. In other words, we do replace at least 2B boundary condition 
points near the edges with the classical ones. Although it may happen that a few yf pile up, but 
this occurs away from the edges. The key property of the family yf is the following bound 

#{« : yf e J} < CN\J\ (5.4) 

for any interval J such that \ J\ > cN^^ and c\J\ < dist(J, [yL,yL+K+i]) < \ J\/c with some small 
constant c. 

Consider the regularized local relaxation measure, which is defined as the probability measure 

/.f ^(dx) = Z-ie-^«""dx (5.5) 
of K ordered points x = {xl+i, . • • , xl+k) in yL+K+i], with Hamiltonian 

^(x) := f - 1 E i°gi^^- - ^.1 +5]Qr(xo, (5.6) 

iG7 i.iei iel 

with a quadratic confinement Ql{x) = (2T(y))^^(a; — 0^)^ as in (4.34) and T(y) = f/s(y)^. The 
potential is given by 

V;n{x) = V{x)-^Y.^og\x-yf\-^ ^ \og\x-y,\ for L + l<t<L + 4B 

j<L j>L+K+l 
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V;s{x) = V{x)-^Y.^og\x^y,\-^ ^ \og\x-y,\ for L + 4B + 1 < z < L + K-AB 

j<L j>L+K+l 

and 

y;«(x) = y(x)-|^^log|a;-y,|-^ ^ \og\x-yf\ for L + if-4i3 + l< z < L + A'. 

j"<L j>L+A'+l 

In other words, we replace the boundary condition y with y^ for the points Xi with L + 1 < i < 
L + AB at the lower edge and similarly for the other edge. The boundary conditions for the middle 
points Xi with L + AB + \<i<L + K — AB remain unchanged. Recalling (4.5), we have in 
particular 

^(x) -H;(x) = § E E [- 1^''^ - Vf\ + log - V^] + (Upper edge), (5.7) 

L-iB<j<L L<i<L+4B 

where (Upper edge) refers to an analogous term collecting interactions near the upper edge. 

Lemma 5.1 Let y G Q . The relative entropies of the measures fiy and /^y satisfy 

S{fil\ti^n + S{^i^'^\^il) < CB''\ogN. (5.8) 

Proof. We start with the following lemma that estimates the relative entropy of any two 
measures: 

Lemma 5.2 Suppose ^i{dx) = Z~^e^^^dx, i = 1,2 are probability measures with Hamiltonians 
Hi on a common measure space. Then 

Si^lM < K^.dH2 - Hi]+E,,,[H, - H2]. (5.9) 

We also have the inequality 

E^, [H2 - Hi] < log Zi - log Z2 < E^, [H2 - Hi]. (5.10) 
Proof. By Jensen inequality, we have 

< SiiiM = j d/ii log (^^^ = J dMi[i?2 - Hi] + log J' 

/dr 

<E^,[H2 - Hi] +E^,[Hi - H2]. 
This completes the proof of Lemma 5.2. q 
Hence we have S{fiy]fiy''^) < I3ili, where 

Qi := (E^r - E^B.r^ E E - vf^ + - y^)\ + (Upper edge) . 

L-iB<j<L L<i<L+iB 

(5.11) 
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Using that Xi ~ yj ^ 1, we clearly have 

C!i < - J2 J2 ^ogix, - yf) + Wy-^ log(x, - y^)] + (Upper edge) 

L-iB<j<L L<i<L+iB 

< CB'^ logiV - B'^E^B.. log(.TL+i - yi) + (Upper edge). (5.12) 



In the first term wc used the trivial estimate Xi — y^ > — — ciV^^ for any j < L. The 
second term will be estimated by Lemma 5.3 below and this completes the estimate for S{Hy\^y''^). 
The other relative entropy, S'(/iy '"^j/iy) can be treated similarly and this proves Lemma 5.1. q 

Lemma 5.3 Suppose r > N~'^ , then for any p > 1 we have 

E^j|log(xi+i < CplogTV (5.13) 

and the same estimate holds w.r.t the measure fJ-y '^ ■ 

Proof. We will need that 

P^j(.Ti+i -yL< N-\) < Cr (5.14) 

for any r G (0, 1). Then (5.13) follows from integrating in r from to 1 and treating the regime 
XL+i -yL> trivially by using xl+i - yL < Vl+k+i ~Vl< CK/N < 1. 

The estimate (5.14) can be proven essentially in the same way as (4.11), just the potential 
^V{xj) of the j-th point in that proof is replaced with ^V{xj) + Q'j{xj). The final estimate 
is somewhat weaker since now the bound on the constant Cy defined in (4.12) deteriorates to 
Cy < Ct~^ < CN. This accounts for the change from N^"^ to N^^ in (5.14). The argument for 
the measure jiy''^ is analogous and this proves Lemma 5.3. q 

5.2 Regularization does not change spacing statistics 

Given that the local relaxation measure and its regularized version /iy are close in relative 
entropy sense, the next proposition shows that their local spacing statistics coincide. 

Proposition 5.4 Let y G G , t = T(y) = T/s(y)^ and assume that for the parameters B = , 
K N'' and f = N^* it holds that 

l + 2b-t-k<0. (5.15) 

Then 



(5.16) 



tei 

as N ^ oo for any smooth and compactly supported test function G. 
Proof. Since Lemma 5.1 and (5.15) guarantee that 

ivs(,M .p. ^ CiVgV ^ ^ 

K K 

with some e' > 0, Proposition 5.4 is a direct consequence of the following comparison lemma which 
was first stated in a remark after Lemma 3.4 in [14], sec also Lemma 4.4 in [15]. q 
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Lemma 5.5 Let G : M — > M &e a bounded smooth function with compact support and let a sequence 
Ei be fixed. Let L be an interval of indices with \I\ = L\ . Consider a measure ui with relaxation time 
T and let qdu) be another probability measure. Then for any ei > Q and for any smooth compactly 
supported function we have 



and 



J^JlJ G{Nix,-x,+,))[q^l]du 



< C 



N^+^^S^iq)T 
K 



(5.18) 



(5.19) 



where Si^{q) := S{quj \ uj). 

Proof. Let q evolve by the dynamics dtqt = Cqt, where C is the generator defined by 

j -fCfdu; = DM) = 2^ / |V/Pdw. (5.20) 

Let Ti = N^'^T. Since q^--^ is already subexponentially close to uj in entropy sense, S'ij((j'ri) l£ 
C cxp(—cN^^)Sa;{q)^ and the total variation norm can be estimated by the relative entropy, we 
only have to compare q with q^^ . 

By differentiation, we have (the summation over i always runs i € I) 

/ ;|^G(7V(x,;-S,))g,,dt^-y" ^ G(iV(a;, - S,))gdcj (5.21) 

i i 

= J^' J j^Yl ^«G'(iV(x, - £;,))) d.qsduj. (5.22) 

Here we used the definition of C from (5.20) and note that the 1/iV factor present in (5.20) cancels 
the factor N from the argument of G. From the Schwarz inequality and dq = 2y/qdy/q, the last 
term is bounded by 



§^l"dsJj2[d.G{N{x,~E,)) 



2 


1/2 




(7s do; 











1 



1/2 



< c 



NSUq)ri 
K 



(5.23) 



by integrating dsS^(qs) = —AD^^{,Jqp). This proves (5.18) and the proof of (5.19) is analogous, q 
5.3 Accuracy of block averages 

In the next Section 5.5 we will compare the regularized local relaxation measure My with the 



reference measure cr^ ' in Dirichlet form sense. As a preparation for this step, we give an estimate 

\B] 



'e 

on the location of the block averages x^^' . Recall their definition 

? + l ^ 



2B 



Xk 



\k-j\<B 
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for any j E + B + 1, L + K — B^. The following lemma shows concentration on a scale C, for 

x^^' w.r.t. /iy and fJ-y'^ ■ The scale C is larger than S but will be smaller than K/N , the length of 
configuration space interval. Thus that the accuracy of the position of Xj decreases from 5 to C, 
but the accuracy of yk is still 5. 

Lemma 5.6 Set ( = t = 2d-l-ip and fix y eQ. For any j e fL + B + 1, L + K - Bj we 

have 



'4\ 



and 



provided 



> C) < cie-^^^' 
> C) < cie-^^^' 
b ip , k 



z < — 03 + mm id , d 1 — 

V 2 2' 22/ 

/or some e' = s'{d, ip) > depending only on d and ip. Furthermore, we have 



nr [B] [B] 



<5C, 



Kj.tX 



[B] 



■ 7 ■ 



<5C, 



7. 



[B] 



<5C. 



Proof. We will need two standard inequalities from probability theory. The first one is 

1 



(5.24) 
(5.25) 

(5.26) 
(5.27) 



M) ■ log 



< log2 + 5(Ai|z^) 



P.(A) 

for any set A and probability measures fx, v. This can be obtained from the entropy inequality 

/d/i < S{ii\v) + \og j e^Av 

by choosing f{x) = b ■ 1^(2:^) with 6 = — logP''(^). Using Lemma 5.1 we thus obtain 

log2 + CS2 logiV 



(5.28) 



^ - -logP^.(A) 



(5.29) 



The second inequality is a concentration estimate. Suppose that the probability measure uj 
satisfies the logarithmic Sobolev inequality (LSI), i.e. 



(5.30) 



holds for any / > with j fduj = 1. Then for any random variable X with E^^X ~ and any 
number T > we have 



E^e^-^ < Ec^ cxp 



(5.31) 



Since the Hamiltonian Hy is convex with V^T-Ly > r ^, by the Bakry-Emery criterion the 
measure /Xy ~ exp{—N'Hy) satisfies (5.30) with Sobolev constant Cs = 2t/N, i.e. 



S{iy I < ATD{y I ify) 



(5.32) 



for any probability measure v (recall that the definition of the Dirichlet form (4.17) contains a 
1/27V prefactor). The same statements hold for the regularized measure Hy'^ ■ 
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For L + B + l<j<L + K-B define the event 

A^Aj = {\ xf^ - E,,j I > C} , with C^N- 



(5.33) 



with a parameter z € (0, 1) chosen later. Using (5.31) for X = ^{^f ^ Epxa;^) and noticing tliat 
|VX|2 = (2B + l)-\ wc obtain 

F^.(A)<2e-5^^^C=r-^ 



Using now (5.29), we get 
assuming 

Using i = 2(i — 1 — we need 



CBt 

P B,r{A) < — 

fly- y I - 



b-t + 2z~l<0. 
b if 



z < d 



Under this condition we have from (5.34) that 



V,.[^\xf-E,.xp\>C)<2e 
Since the measure /iy is also concentrated by the LSI, we have 
P B.r (\xf^ - E n,.xf^ \>C) < 2e-^' ^ 
and together with (5.35) we have 



2 2 



IB]\ 



(5.34) 
(5.35) 
(5.36) 
(5.37) 

(5.38) 



EB..xf~E^.xf 



(5.39) 



r D] 

Therefore X; is concentrated on a scale Q around the same point w.r.t both measures /iy and 



Using (5.32) and that y G we get 



(5.40) 



Hence by (5.28) and (5.34) we obtain 



■ My 



(A)<i^^^^±i=i=<^^o 

-log (A) ^ Be 



provided that 



, k b 

z < d h -. 

2 2 



(5.41) 



(5.42) 
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Now by the definition of y e in (4.7) we have 



(l 



SB] 



[B]\ 



< 



< 



I V E 

(2B + 1)C ^ 



{B\\1 



1 



E 



(2B + 1)C^ ^ '"^^ 



I |2 <- 



2 



using (5.37). Combining it with (5.41) we obtain 



E^y 



{B\ 



<2C. 



(5.43) 



Finally, since y G tj, we have 



lE..y4 



which, combined with (5.39) and (5.43), yields (5.27). This completes the proof of Lemma 5.6. q 
5.4 Proof of Proposition 4.5 

Now wc will compare the regularized local relaxation measure with the reference measure 

CTg'* in Dirichlet form sense. Recall their definitions from (5.5) and (4.33), respectively, and recall 
that T = T(y) f/.s(y)^. Here s = s(y) is a function that is approximately 1 for good external 
configurations y £ G (sec (4.23)). 

The result is the following comparison of local gap statistics. Combining this result with 
Proposition 5.4 and checking that the condition (5.15) is satisfied with the choice of parameters 
given below, wc arrive at the proof of Proposition 4.5. q 

Proposition 5.7 Fix ip < j^. Let y E G ■ t = T(y) = f/s{y)'^ and assume that for the parameters 
S = N-"^, B = N^, K = N^ with d=l-(p,b^8ip, fip. Then with t = 2d-l-ip=l-3(p 
let T = N^* with t := 2d — I — ip = I — 'iip. Then 







(5.44) 



as —> oo for any smooth and compactly supported test function G. 



Proof. The key technical estimate is the following lemma whose proof will take up most of 
this section. 

Lemma 5.8 Let <^ > 0. Suppose B = , K = with < b < k < 1, and S = N^'^ with 
d G (0, 1). Suppose that these parameters satisfy 



. b ip k b\ 

1-b < -ip + mm[d - -- + 



(5.45) 



i.e. one can choose a number z > 1 — 6 and satisfying (5.26). Let y Cz G ~ Gs,£o good 
configuration. Assume that £o !i where e' = s'{d, ip) is obtained in Lemma 5.6. Assume that 
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the equilibrium measure pv is away from the edges. Let t = N ^ with t = 2d — 1 — (p. Then 
the Dirichlet form of Py '^ with respect to the reference measure is bounded by 



T 

K 



I 4'^) < Cf(logiV) + + ^ + + iV^/^+^J + Cle--^^'^^ (5.46) 

The prefactor t jK is for convenience; the local gap statistics of two measures are approximately 
the same if tD/K 0. More precisely, we have the following general theorem which is a slight 
modification of Lemma 3.4 [14] (see also Theorem 4.3 in [15]). This result was originally proven 
for /3 > 1, but by a regularization argument it extends to any /3 > 0, see Lemma A. 2 of [17] for 
details. 

Lemma 5.9 Let G : M — > K. &e a bounded smooth function with compact support. Consider a 
measure on Yik {x : xi < . . . < xk} C defined by 

dc^ - e-'^^^dx, H(x)=Ho(x)-^ ^og{xj-x,), (5.47) 

l<i<i<K 

with the property that W'^T-Lq > holds for some positive constant t. Let qdoj be another prob- 
ability measure. Let I C {1, 2, . . . , K -~ 1} be an interval of indices. Then for any ei > and for 
any smooth compactly supported function we have 



I^E/ G(A^(^. - [g - l]dc.| < C^ ^'^^l^^^ + Ce-^" VKM, (5.48) 

where D^{y/q) D(qu) \ w). 

We will apply Lemma 5.9 for the measure w = cxg'*. It has the form (5.47) except that (Tg''' 
is restricted to the interval [yLTVL+K+i], i-e. the relations j/l < xl+i and x^+k < Vl+k+i also 
hold in addition to the ordering relation xl+i < xl+2 < . . . < -^l+a'- Notice that the Hamiltonian 
"Hg'^ of the measure ag'^ contains a term [log(a;L_|-i — ?/l) + \og{yL+K+i — xl+k)] ■ This term 
confines the particles in the interval [yL, Vl+k+i] exactly as the term log(xi+i — Xi) guarantees 
the ordering constraint Xi < Xi^i. Hence the regularization argument in Lemma A. 2 of [17] can 
be used to treat the additional constraints, yL < xl+i and xl+k < Vl+k+i- 

The proof of Proposition 5.7 now follows from Lemma 5.8 and Lemma 5.9 with lu = tTg'* and 
qduj = l^y'^ ■ The parameters b,k,d G (0, 1) have to satisfy the following relations from (5.45) and 
from the requirement that the right side of (5.46) converges to zero: 







b 


< 


k 






1 


-b- 




< 


d- 


b 

2 ^ 


2 


1 


-b- 


f <^ 


< 


d- 


k 

2^ 


b 
2 


1 - 2d + (/3 4 


-2k 


- 1 


< 









l-2d+{2d- 


- 1 - 


^) 


< 









4fc - 3 + (2d - 


- 1 - 




< 









i + ip-2d + 3 


-b 


- k 


< 









l-2d + 


2lp- 


3 

^5 


< 


0. 







39 



It is easy to check that aU these conditions are satisfied if, e.g. 

39 1 

d=l-Lp, h = %ip, k = —ip, < < — . 

This choice is not optimal for the above system of inequalities, but we took into account that 
the parameters will also have to satisfy (5.15) so that we could combine Proposition 5.7 and 
Proposition 5.4 to arrive at Proposition 4.5. 

Finally, the entropy term ^(/iy ''^ | crl'^) in (5.48) can be estimated by the Dirichlet form via 
the logarithmic Sobolev inequality. This completes the proof of Proposition 5.7. r-i 



5.5 Dirichlet form estimate: proof of Lemma 5.8 

By definition, 

^ e L+l<j<L+K 

where Zj is defined as follows: For L + l<j<L + AB, we set 



k<L-2B J fc<L-2B J K 

k>L+K k>L+K 

(recall that 9'^ = 9j/s and we set Ws{x) = s^x^). Note that the summation at the lower edge is only 
for fc < L — 2B instead oi k < L because the interaction terms near the boundary cancel by (5.3). 
Moreover, notice that the linear terms, coming from the derivative of the quadratic confinements 
(see (4.3) and (4.34)), cancel each other 



T ■' T T 

by the choice of T(y) = f/s(y)^. 

Similarly, for L + K - AB < j < L + K , we set 



k>L + K-\-2B 3 
k<L 



Finally, for L^4:B<j<L-\-K- AB, we define 



k<L J k<L J 



k>L-\-K + l k>L + K + l 



Notice that here yk is not replaced with since only interactions for xj 's near the edges have been 
regularized. Moreover, the interactions with the boundary terms yk, with k = L and k — L + K + 1 
cancel out since yL — and j/l+k+i = ^l+k+i matching construction. 



Now we estimate the size of Zj in each case. 
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Case 1: L + AB < j < L + K — AB. The first step is to decompose Zj as 

5 



(5.49) 



where 



n] 



pv{y) 



Xj - y 



pw, (y) 



Xj~y 



^2 ^ ^2,(0,.- ^ ^2,up 



^ 6^ - vk 



3^low ^ ^3, up 



pvjy) 

Xj - y 



dy 



N ^ 



1 



k>L+K+l 



pviy) 



^ Xj yk Jy^+K+i ^3 y 



dy 



Vl+k+1 



PWs iy) 
00 -^j ~ y 



dy 



- y ^— 

N ^ X,- 9', 

\ k>L+K+l J « 



pws jy) 

Xj - y 



dy 



pvjy) ~ pwAy) 

Xj-y 



dy 



(5.50) 



Here we also used that [1/^,2/^+^+1] = [^i,, when establishing the limits of integrations. 

By the equilibrium relation (3.2) between V and pv, we have 



From (4.25), we have 



P^f = C 



(n,-e'j? ^ c 



(5 H r 



(5.51) 



(5.52) 



Since pv € away from the edge, and so is the semicircle density pw^ 1 we have by Taylor 
expansion 



VL+K+l 



Pviy) - pwAy) 



dy 



< 



VL ^3 y 

pvjxj) - pw, {x,) + 0{xj - y) 

VL 



(5.53) 



Xj -y 

<C[\ \og{x.j - yL)\ + I log(yL+A'+l - Xj 



dy 



K SN 



Here we used (4.27) and (4.30) and the fact that pw^ (x) — pwix) = 0{\s — 1|) away from the edge 
together with (4.23) to estimate 



\pv{x) - PWs{x)\ < C 



K SN 
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for any x € [yLTllL+K+i]- The logarithmic terms after taking square and expectation w.r.t. wiU 
give rise to an irrelevant logiV factor by using Lemma 5.3 

E^B,. [I log(x, - vl)\ + I \og{yL+K+i - < Clog TV. 



We now estimate the main error and we will deal with the first term only, coming from the 
lower edge, the second one can be treated similarly. We write it as 

pv{v) 



k<L 



-oo ^'j y 



-dy 



2,3 



with 



n 



2,2 
3 

2.3 



Ay 



Pv{y) 



Pv{y) 



dy 



71, ^3 y 

N ^ La;,- - 



- Ik Xj ~ yk 



With C, = N ^ with z is given in Lemma 5.8, define the event 

A= -7fl| <6C, lL + i3 + l,L + A--Bl}, 

then its complement has very small probability, 

from (5.24) and (5.27). On the event we simply estimate 



IN 



k<L J -oo Xj y 



< 



yL - VL-i 



+ C log(a:j -yi) 



therefore 



Pv{y)_^y - 



(5.54) 



< 



+ IE.-- I ^og{x, ~ yL)\') '^'(P,B,. {A^)Y^' < cie--^^'^^ (5.55) 



by using Lemma 5.3 and — J/L-il > exp(— A^^°) from y £ G. Here we used that Eq < e'/lO. 

Now we continue the estimate on the set A and we consider the three terms in (5.54) separately. 
For the first term we write 



^fc<L^^^-^^ 4 ^3 



Npviy) 



dy 



1 

N 



E 



T'^Npviy)dy 
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where we have used that J^^^^ Npv — 1 and 

r 



Ik 

k ik-y 



" {xj-y)(xj--iky 



Recall that L > kN 3> 5 and Xj e [y^, j/l+k+i] = [iltIl+k+i] + 0{5). For /c < ^k-/V we know 
that |7fc — I > c with some positive constant. Hence we have 



1 r< pk+i n 



k<KN/2 /£<re7V/2 k<KN/2 

since 7fe+i — 7fe < CN~^/^k~^/'^ near a square root singularity of py at the edge. For the regime 
k > \kN we can use |7fc+i — 7fc| < CN~^ to get 

_ \ < — > 

N ^ ^-N ^ N(xi--ikY 

KN/2<k<L KN/2<k<L ' ■' ' ' 

C 

< 



- E 



KN/2<k<L y-^j-B Ik) 

C 1 c 

Here in the second inequality we used that on the set A we have 

Xj > xfJs > l\% - 6C > 7j-2i3 - 6C > 7L + cBN-^ > jk + cBN-^ (5.56) 

ioT k < L using j > L + AB and thus 7j-2B — 7l > cBN~^ ^ 6C, since z > 1 — h. Therefore 
Xj — 7fc > x -_^ — 7fe > 0. In the third inequality we performed the summation and used that 7/c 
is regularly spaced. In the last inequality we again used (5.56). In summary, we have shown that 

on the set A and we have seen that the contribution from A'^ is subexponentially small (5.55). 



Now we consider fi^'^ on A. We have 



dy ^ C8 



\^Y\ < C / < _ , (5.58) 



by using Xj-^l > lj-2B-lL-K > c{^j-2B-"/L) from (5.56) and from7j_2B-7L > cBN~^ > 6C, 
moreover xj ~ yL > xj — — S > c(7j_2B — 7l) by \-fL — 2/l| < ^ (from y G ^) and S <C BN~^ 
(from (5.45)). Thus 



B 

L+4B<j<L+K-4B 
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For the third term ft^'^ we have 



L+iB<j<L+K-iB 



(5.59) 



-fly ' 



E 



< 



L+'iB<j<L+K-iB 

■1(A) E 



-E 



k<L 



1 



L+AB<j<L+K-iB L fc<L 



-y- 



Xj - Vk Xj - 7fe 

H 2 



(yfc - 7fe) 



{xj - yk){xj - jk) 



We spht the summation over k into two terms: kN/2 < k < L and < kN/2 and separate by a 
Schwarz inequahty. 

First we consider the case kN/2 < k < L. Expanding the square, we need to bound 



E B.^l(A)^ 



E E 



E 



\Vk -IkWVa -7q| 



KN/2<k<L KN/2<a<L L+iB< j<L+K 



{Xj - yk){Xj - lk){Xj - ya){Xj - 7a) 
1 ^r^ 1 



(5.60) 



E 



KN/2<k<L 



L+4B<j<L+K ^""^ y'^)^ KW/2<a<L (^^' '^"^^ ' 



where we used another Schwarz inequahty and the factor 2 accounts for a similar term with the 
role of k and a interchanged. 

In the case kN/2 < k < L wc have [7^ — yk\ < S. Then (5.60) is bounded by 



2^2 



E 



E 



■1(A) E 



1 



E 



1 



L+4B<j<L+A' fe<L 

<C6^ 



i^j -ykf (^j ^ 



(5.61) 



E 

L+iB<i<L+K 



E fl,xl(A)- 



< C5' 



Here we used 



-y 

a<L 



1 



^Xj 7a) 2 



< 



C 



< 



c 



Xj ~ 7L Xj 



VL 



relying on the regularity of 7a and using, from (5.56), thatxj— 7a > a;j— 7l ^ cBN~^ which is much 
larger than the spacing of order iV~^ of the 7-sequence. In the last estimate Xj — jl ^ \jl ^ yL\ 
was used (since BN~^ ^ 5). Similarly we could perform the k summation 



E 

k<L 



< 



c 



(xj ~ykf Xj - yL 



since Xj - yk > xj - •jk ~ S > c{xj - 7fe). 

To perform the j summation in (5.61), we use 



1 



< 



and then we recall that apart from a set of subexponentially small probability, we have 
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from Lemma 5.6. Since C ^ BN ^ and x^^}b " Vl > cBN ^ from (5.56), we see that 
V 1 ^ C ^ CN ^ CN^ 

L+4B<j<L+K ^ ^ L+iB<] \lj-B VL) IL+3B VL 

On the exceptional set one can just use the trivial bound {xj — j/l)^^ < C!{xj — 11)^"^ < CN^B~^ 
from (5.56). 



Consider now the case k < kN/2 in (5.59). We have 

E 



(\-)_}_ "V "V 'V Ij/fc - IkWya " 7a| (^P,9\ 

^ I (x, ~yk)(x, ~jk){x, ~ya){x, -Ja) ^ 

k<KN/2 a<KN/2 L+iB<j<L+K ^ ^"""^ '"^^ ^"-'"^ ^ '°-> 



k<KN/2 L+iB<j<L+K ^ ^ a<KN/2 ^ ^ ' ' 

<— E^.,.1(A) E k-7feP<CA'iV-2/5+^, 



k<KN/2 



where we used that all denominators are separated away from zero and Lemma 3.6. Furthermore, 
in the last inequality, we have used Lemma 3.6 for k > N^/^^"^ and we used \yk — 7fc| < 0(1) for 
k < N^/^+'^ from y ^ Q and (4.7). Similar comment applies to all edge terms in this proof and we 
will not repeat it. 

Summarizing, we have shown that 

E^B,.1(A) E [nYY <^^^^ + CKN-^/''+^. (5.63) 

L+iB<j<L+K-iB 

Finally, we need to estimate ilj^ in (5.50). It can be treated exactly as 57^'^ and the result is 

on the set A and the contribution from A'^ is subexponentially small as in (5.55). 

Case 2: L < j < L + 4B. (There is a third case L + K - AB < j < L + K which is identical to 
Case 2 and will not be treated separately) . We decompose Zj as before and the only modifications 
are 

J ' N ^ Xi - 6' x.~y ^ 

k<L-2B ■' k -J oa 1 y 

^4 r^''^' pv{y) - PwAy) ^ f"'-'^ P^Mdy 

'' JyL-2B ~ V JyL-2B ~ V 
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We now estimate the main error term r2^''°™, and we write it, as before 

^2,lon, ^ ^2,1 ^2,2 ^2,3 

7 7 ' 7 ' 7 



with 



,2,2 



-{- E — 

■dy 



k<L-2B ■> ' 

■y---^ pv{y) 



pvjy) 

Xj-y 



dy 



'iL--2B ^3 y 



,2,3 1 



- y \^ 

N ^ Vx, - 



k<L-2B 



Xj - 7fc Xj - yk 



We have 



- E 



y-ik 



fe<L-2B-'Tfc 



{xj - y){xj -7fc) 



Npv{y)dy 



< 



C 



- y 



7fe+i - 7fc 



k<L-2B ^ ■> 
C C 



{xj - 7fc)^ 

C C 

|7fc+i-7fc|< 

k<KN/2 

using that Xj > yL Jl S > jl~2B + cBN^^. The estimate of il^'^ is trivial 
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(5.65) 



KN/2<a<L-2B 



{Xj - 7a) 2 



TV 



k<KN/2 



where we again split the summation over k into kN/2 < k < L — 2B and k < kN/2, yielding the 
two terms, similarly to (5.61) and (5.62). 

The estimate fij^''""' is analogous to that of ^^'^ ■ The first term of fi^ is estimated as before in 
(5.53). The additional second term in fJ^ is trivial by recalling |7l-2S — ^l-2bI — CB'^N~^ + C6 
from (4.25): 



' VL-IB 



Pw. jy) 

Xj-y 



■dy 



< 



yL-2B - 0'l-2b\ ^ C6 + CB'^N-^ ^ CSN CB 



cBN-^ 



cBN-^ 



< 

- B 



N ' 
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since the denominator can be estimated by using Xj — yL-2B > Vl — 2/L-2B > 7l ~ 1l-2B — 2(5 > 
cBN^^ and 

- 0'l-2B >yL- 0'l-2B >1L- 1L-2B - 25 + (7i_2i3 - 0'l^2b) > cBN'^ 

where we used 6 < BN^^ and B < A^. 

Collecting all the error terms into (5.46) and removing some redundant terms, we have thus 
proved Lemma 5.8. q 

Acknowledgement. We would like to thank M. Ledoux for pointing out an error in the statement 
of Lemma 3.9 in an early version of this paper. 
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